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Abstract 

We consider a 2D electron system in a strong magnetic field, where the local 
Hall resistivity Pxyi^) is a function of position and Pxxi^) is small compared 
to Pxy Particularly if the correlations fall off slowly with distance, or if fluc- 
tuations exist on several length scales, one finds that the macroscopic longi- 
tudinal resistivity Rxx is only weakly dependent on Pxx Q-nti is cipproximatGly 
proportional to the magnitude of fluctuations in pxy This may provide an ex- 
planation of the empirical law Rxx oc B where Rxy is the Hall resistance, 

and B is the magnetic field. 
PACS: 73.40.Hm,73.40.-c,64.70.Ak 
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Almost ten years ago, it was noticed that the longitudinal resistance Rxx of a quantum 
Hall system looks very much like the derivative of the Hall resistance R^y with respect to 
filling fraction [jT|. More recent experiments [0] have shown that in high mobility samples, 
the relation 

Rxx = (1) 

with B the magnetic field, and a a constant, holds amazingly well over a wide range of tem- 
peratures and magnetic fields, including both the integer and fractional quantized regimes, 
and unquantized regions where Rxy varies linearly with B, as for a classical system. The 
constant a is sample dependent, and varies somewhat with temperature, but is typically in 
the range 10^^ to 10^^. Despite the simplicity and generality of this empirical relation, it 
has defied explanation for almost a decade. In this paper, we propose an explanation based 
on inhomogeneities in density over long length scales. 

We consider a two dimensional electron system where the local electron density n{r) is 
fixed by the charged impurity distribution in a nearby doping layer [Q. We assume that the 
fluctuations in the density 6n{r) = n{r) — hq are much smaller than the average density hq, 
and we define a dimensionless correlation function g such that {6n{r)6n{r')) = fj^'^g{\r — r'|) 
with g{0) = 1, and g{x) ^ as x — > oo. Although we will later be concerned with the effects 
of disorder over several length scales simultaneously, to begin with, we will consider the case 
where there is only disorder on a single length scale a. As an example of a correlation 
function with a single length scale, we can consider g{x) = exp(— [^]^). Typically, the length 
scale a of any disorder will be greater than or equal to the setback distance of the doping 
layer. If a is sufficiently long, we can define a local resistivity tensor P (r) which depends 
on the local conditions. We assume that at a fixed temperature T, the local Hall resistivity 
Pxj/(r) = —pyx{r) is purely a function of the local filling fraction z/(r) = (j)Qn{r)/B with 
00 = hc/e. This assumption is certainly valid at high temperatures, where pxy takes on 
the classical value, Pxy(r) = e^z/(r)//i. The assumption should also be valid at least as 
a first approximation in the ranges of T and B where the integer or fractional quantized 
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Hall plateaus begin to develop. Then, if the density fluctuations are not too large, we may 
write Pxy(r) = f{^{r)) = Po + 5p-^y{r) where po = /(z/q), z/q = (poUo/B, and 5p.j:y{r) = 
f'{^o)4>oSn{r)/B. Thus we have {6p^y) = and {6pr,y{r)6pr,y{r')) = X^gUr - r'|), with 
A = f'{vo)(f)op/ B. We also assume that there is a local diagonal resistivity pxx(r) which may 
depend on the temperature and the magnetic fleld, but which is small compared to A. This 
assumption should be quite reasonable for high mobility samples. For simplicity, we shall 
take Pxxij) = P, independent of r. 

We now want to flnd the macroscopic resistivity Rxx for the system. The current proflle 
j(r) is defined by Maxwell's equation V x E = and current conservation. Using E = p j, 
we can thus define the current with the system of equations 

Vx(pj) = , V-j = 0, (2) 

along with appropriate boundary conditions and the condition that the spatial average of 
j have a specified value jo. By using the second equation of (0), we can rewrite the first 
equation as 

(Vpxy X j) + p(V X j) = 0, (3) 

which is clearly independent of po. The total dissipation per unit area is given by ^-RixljoP- 
Since the local dissipation |p|jp is independent of po, Rxx must also be independent of po. 
On dimensional grounds (since Rxx depends only on A, p, and the function g), we expect 
that in the limit p — » 0, the longitudinal resistance should have a scaling form 

Rxx = CA^-'^p- (4) 

where C is a nonuniversal constant that depends on the form of the function g. Calculations 
described below suggest that the exponent uj can be close to zero, particularly if we consider 
the effects of disorder on several length scales simultaneously. If indeed we find that u; ~ 0, 
then we have Rxx ~ CA = C f'{vQ)(f)op/ B. On the other hand, the Hall resistance will be 
given by Rxy = Po + ^(A^) ^ /(z/q) so B^ = f{uo)<Pono/ B, thus establishing Eq. (|Tp 
with a = Cp/riQ. 
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A different model, which may be more apphcable for high mobihty systems at temper- 
atures that are not too low, has dissipation arising from electron-electron scattering, which 
causes no dissipation for a uniform flow velocity. If the magnetic field is not too strong, we 
may represent this by an electron fluid viscosity rj, instead of a resistivity p, so that the local 
dissipation is given by T]n'^\V x rather than with j the local current density. In 

this case we find the similar result Rxx = Ca^'^^X^^'^r]^ . If the exponent u is sufficiently 
close to zero, so that we may neglect the field-dependence of 77, we may again derive the 
resistivity law (j^). 

The conductivity problem described above has been studied by various authors 0,^,0. 
These analyses conclude that the exponent u for the case of a local resistivity p, with finite 
range correlations in p^y, is = ^ ~ 0.23. A similar analysis for the case of a local viscosity 
gives = ^, as will be discussed below. (This result has also been derived by Fogler and 
Shklovskii [0.) 

A starting point for understanding this conductivity problem is the observation that in 
the limit p — 0, the current j(r) is confined to a contour of constant Pxyi^^) which percolates 
across the system [^^. If 6pxy{r) is statistically symmetric about zero, then this percolating 
contour |g] is defined by Pxyij) = Po- As p is increased, the current spreads out somewhat 
from this contour. The physical interpretation of u; ~ is that (in either the viscous case 
or the resistive case) as we decrease the dissipative part of the local resistivity (fixing the 
potential difference across the system), the current profile readjusts such that the total 
dissipation remains approximately fixed. 

If we imagine one of the percolating contours of Pxyi^^) = Po to be locally pointing in 
the y direction, we can write locally 6pxy{x,y) = xQ{y) for some function Q. Since we 
suspect that the current will follow along this contour, we propose that the y component 
of the current can be written in the form ]y{x,y) = JS[x/w{y)]/w{y) where J is the total 
current carried by the channel, w{y) is the local width of channel, and 5* is an unknown 
function satisfying the normalization condition J dzS[z] = 1. Of course, once we know j^, jx 
is determined by current conservation. In the small w limit it is easy to show that a solution 
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to the system of equations (|^) is given for the functional form S[z] = exp(— 2:^/2) / V^tt when 
the width w of the channel is determined by the differential equation 

d[{wQr]/dy = -ppQP-' (5) 

where p = 2 here. In the case of viscous dissipation, we find a differential equation similar 
to (H), except that p is replaced by r], and p = 4. In the viscous case, the function S must 
be chosen to satisfy the differential equation d^S{z)/dz^ = —zS{z). 

As an explicit example, we consider the periodic geometry 

TT{x+y) \ ( n{x-y) 



5pxy{r) = A cos ^^ll£±yl^ - cos 



which has a square lattice of saddle points at 
(^, |) G Z^, and we take p ^ A ^ po- We assume an average current jo in the y direction so 
that the current flows down channels of fixed Pxy{T^) = Po along the x = integer lines. Near 
the channel at x = 0, for example, we write 6pxy{x,y) = xQ{y) with Q{y) = —Hsin^iry/s), 
where H is the characteristic perpendicular slope of the function 6pxy In this geometry we 
have H = 27iX/s, but in a more general geometry, H ~ A/a, where a is the characteristic 
length scale of the density fluctuation. 

Integrating w from one saddle to the next in Eq. (^), we find that the width of the 
channel is given by w{y) = ps/ttH\ sec(7r?//2s)| in the resistive case. Of course the small w 
approximation breaks down near the points where this expression diverges, so the solution 
must take a different form very close to these points. The important thing to extract 
from this result is that the characteristic width W of the channel scales as ly ~ ps/ H. 
Numerical solutions of the system of equations (H) to find the exact current profile in 
this periodic geometry (as well as in other simple geometries) support this conclusion. Since 
we have determined the current profile, we can also calculate the local dissipation p|jp and 
thus extract the resistance of this current carrying channel. In general, we find that the 
resistance of a channel of length s scales as Rs ~ ps/W ~ HW . A similar argument for the 
case of viscous dissipation yields the results W ~ {r]s/H)^ and Rg ~ r]s/W^ ~ HW. 

We must remember, however, that we do not actually have a regular array of saddle 
points. We have instead a tortuous conduction network. We define a percolation problem 



by choosing a cutoff e, considering all points such that \5pxy(j:)\l\ < e to be "conducting" 
points, and all other points to be "insulating." This type of problem has a critical value 
of ec = 0, provided we take 5pxy to be statistically symmetric about zero. Our actual system 
wiU have thin channels of width W centered at the percolating 5pxy = lines, and we will 
require for self consistency that e ~ W/a for these channels. The percolation network has 
a web-like topology consisting of essentially one dimensional conduction paths that snake 
through the system and intersect each other once every distance ^. More properly, the cor- 
relation length ^, is defined as the size of the largest insulating islands in the conducting 
network. As we approach the percolation threshold by letting e — > 0, the conduction chan- 
nels become thinner, connections between conduction paths are broken, and the correlation 
length diverges as ^ ~ ae''^. For the case where the disorder has only one length scale, it is 
known that i/ = | |,|,|. 

Although the conduction path in this percolation problem can branch somewhat 
on scales smaller than ^, forming dead-end loops and small multiply connected regions, 
it is assumed for now that the current effectively follows a one dimensional path around 
the perimeter of these branched region The length of such a path connecting points 

separated by a distance ^ scales as a{^/a)^ where D is the fractal dimension of the path. 
For this particular problem, where there is only one length scale, it is believed that -D = | 

i- 

Due to the web-like topology of the network, we can think of the percolation network as 
an effectively periodic system similar to the one we analyzed above . Here the periodicity 
length is ^. However, since the channels are not straight, the length s of the conduction 
channels must be taken to be the typical path length a{C,/a)^ between intersections of 
macroscopic conduction paths. Thus we have W ~ psjH ~ a^/p^^JXaP . On the other 
hand, our self consistency condition requires W ~ ae. Using ^ ~ ae"'', we find p ~ Xe^^'^^ . 
Since the system is homogeneous on length scales larger than ^, and since there is typically 
only one conduction path per square of length ^ on a side, the macroscopic longitudinal 
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resistivity is given by Rxx ^ Rs ^ HW ~ Ae ~ \^~'^p'^ where u = 2+^0 ~ n '^^^^ 
a single disorder length scale. This result agrees with the previous work of Refs. 0] and . 
An analogous calculation can be performed for the viscous case with a single disorder length 
scale (using W ~ [rfs/H]^ instead of ~ ^Jps/ H) to give the result uj = = ^, in 

agreement with Ref. [0. 

Although these arguments seem quite reasonable, and are in decent agreement with 
existing numerical work 0], it is possible that they are not exact. In particular, we are 
not very certain of the fractal dimension D of the conduction path. For example, if the 
current were to readjusts itself only very slightly so as to cut off dead-end loops, D might 
be changed. However, it is clear that 1 < -D < 2, and this uncertainty allows only a small 
range of values of < < ^ in the resistive case, and ^ < uj < j^'m. the viscous case). 

We now consider the effect of having disorder on multiple length scales. As an extreme 
example, let us consider a system that is disordered on a microscopic length scale a (as 
above) and also disordered on a macroscopic length scale a' such that a <^ a' <^ where L 
is the size of the system. On scales much less than a', the calculation described above should 
hold. I.e., given a local resistivity tensor p such that p^y has some typical fluctuation A, and 
Pxx = P, we can use the above approach to calculate that the resistivity tensor p measured 
on a scale larger than ^ but smaller than a' is given by p'^^ = p' ~ A^p^ with p'^^ given 
by the local mean value of p^y We can then repeat the argument to account for disorder 
on the scale a' by using the tensor p as a local resistivity. Here, we similarly find that 
the macroscopic resistivity is given by R^x ~ ATi(p')A ~ ATS^pifg. Thus the exponent u) is 
effectively squared if disorder exists on two very different length scales. Clearly, this approach 
can be extended analogously to account for disorder on any number of well separated length 
scales, with the result that the exponent lo can become arbitrarily small if the disorder 
extends over a sufficiently wide range of scales - i.e., if the correlations fall off sufficiently 
slowly. Isichenko has previously considered the effect of long range correlated disorder and 
has similarly found that uj can be arbitrarily small for certain correlations that fall off very 
slowly i. 
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We can also show rigorously for a slightly modified model, that the exponent uj can never 
be negative. In the modified model, we take Pi.y(r) = poc^'-'"-', and pxx{'^) = pe^^^\ where /i(r) 
is dimensionless here and has a symmetric distribution about zero. This model is statistically 
self dual, in the sense that for any member A of the ensemble with a local resistivity tensor 
Pa (r), there is another _B, having equal probability, with [p b (r)]^ = [P a ij)Y^{PQ + P^)- 
It can be shown that the ensemble-averaged resistivity tensor R satisfies R=R [pg + p^], 
which implies i?^^ + R^y = (pg + p^). (The proof employs a transformation analogous to 



that used by Dykhne and Ruzin |T^, with the choice a = c = 0, d = b^^ = [p^ + p^j^^^ in 
Eq. [8] of Ref. [[l^.) Therefore, Rxx cannot diverge for p — > 0, and we must have u > 0. 
More generally, one expects that the exponent u will always be > for any model where 
the macroscopic resistivity is isotropic 0]. 

The hypothesis that disorder exists on several length scales seems plausible. Although 
one certainly expects to find disorder on the scale of the dopant setback distance, there is 



also evidence for inhomogeneities in the electron density on the scale of a few percent |]TT 
and the length scale of these inhomogeneities may be relatively large. In this regard, it 
is interesting that the values of axx obtained from the surface acoustic wave experiments 
of Willett et. al. using a conventional theoretical analysis, in the low frequency limit, are 
systematically smaller by a factor of ~ 3 than the values obtained from macroscopic con- 
ductivity measurements on similar samples, over a range of fields and temperatures. (This 
discrepancy was removed by an explicit renormalization of a parameter am in the experi- 
mental analysis |[TT| , p!2[| ). If the sample has sufficiently large density fiuctuations on a length 
scale large compared to the ultrasound wavelength together with smaller fiuctuations on a 
shorter length scale, this could conceivably account for the observed discrepancy. However, 
there remains a problem in accounting for the experimental values of the coefficient a in Eq. 



(0) via large scale inhomogeneities. Analysis by Dykhne and Ruzin |]T0[ in the case where 
Pxyi'^) can only take on two discrete values, pi and p2, gives a maximum value of Rxx equal 
to |pi — P2I/2, when the value of p is small. Assuming this estimate to be an upper bound 
on Rxx for the continuum case, we see that in order to obtain a coefficient a ~ 0.03, as in 
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typical experiments 0, one would need to have large scale density fluctuations of order 6 
percent, which is probably unreasonably large. 

If the exponent uj were precisely zero, then the resistivity law (|1]) would hold independent 
of the microscopic details of the local dissipative resistivity (or viscosity) so long as the local 
dissipation is sufficiently small. If uo is small, but nonzero, then the accuracy of the law 
will be determined to a some extent by the behavior of the local resistivity. In the high 
temperature regime, for example, if the local resistivity p itself varies as with a ^ \, 
then at high temperatures the macroscopic resistivity would vary as B^'^^'^. 

Theoretical results that in various particular circumstances, within the quantum Hall 
regime, the macroscopic longitudinal resistivity is determined by local fluctuations of the 
electron density n{r), or equivalently of the local Hall resistance pxy{v), have been discussed 
previously by a number of authors p,^|10|,|12|,|13[ . Much of this previous work has concen- 
trated on the narrow transition regions between neighboring well-developed Hall plateaus at 
low temperatures. In the present note, we have been more concerned with the regime of high 
temperature, where Pa:j/(r) may be more properly considered to be a continuous function of 
the local fllling factor and hence a continuous function of position in the sample. Moreover, 
we have shown that a percolation analysis of the conductivity problem may be used to ex- 
plain the empirical resistivity law (|l]) which should hold particularly well if disorder occurs 
over a wide range of length scales. Further analysis of the inhomogeneities occuring in ac- 
tual samples is clearly necessary to complete this argument, as is a better understanding of 
dissipation processes at the shortest length scales. 
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